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Abstract 

We consider the sum of dichromatic polynomials over non-separable 
rooted planar maps, an interesting special case of which is the enumer- 
ation of such maps. We present some known results and derive new 
ones. The general problem is equivalent to the g-state Potts model 
randomized over such maps. Like the regular ferromagnetic lattice 
models, it has a first-order transition when q is greater than a critical 
value but is much larger - about 72 instead of 4. 



Introduction 

The dichromatic polynomial of a graph ^ is a polynomial xiG'i^iU) 
in two variables x,y. It has been studied by Whitney ||35[|, Tuttep8[| 
and others. In particular, Tutte obtained a functional relation for 
the generating function <I> of the sum of dichromatic polynomials over 
rooted planar maps. 1 31] 



Fortuin and Kasteleyn||15|, |19|] showed, from the viewpoint of sta- 
tistical physics, that xiG'i^iV) is equivalent to the partition function 
of a g-state Potts model on ty, with g = (x — l)(y — 1). This means that 
$ is the partition function of the sum of maps. From it one should 
be able to obtain the thermodynamic properties of the planar Potts 
model, averaged over all rooted maps. This is not without interest: 
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usually one considers such models on a regular lattice, but physical 
lattices are seldom (if ever) completely regular. They contain defects, 
and one should in fact average over such defects, i.e. over a lattice 
with some degree of randomness. In this sense ^ gives the average 
over over a completely random "lattice" . 

There has also been much work in statistical mechanics and field 
theory on summing over planar Feynman diagrams, using field-theoretic 
techniques and matrix models. 21, |3^, |l|, 39| The 



main difference between that approach and the one used here is that 
our maps are "rooted" , which means that they are given extra weights, 
corresponding to the various ways a given map can be rooted. It ap- 
pears that this is merely a "bounday effect" which does not influence 
the asymptotic large-graph behaviour. 

Here we do try to reduce the randomness by restricting each map 
G to be non-separable (i.e. irreducible). The required analogue of 
Tutte's functional relation for this case has been obtained by Liu (eqn. 
4.17 of 1^^). There are three cases for which the relation can be solved 
explicitly: 



Case 1 
Case 2 
Case 3 



X, y both large, of the same order; 
X, y both small, of the same order; 
q = {x — l)(y — 1) = 1, or X large and y of order one, 
or y large and x of order one. 



We present the solution in all three cases. The first is rather trivial. 
The second appears to be both non-trivial and new: one interesting 
feature we observe is that for non-separable maps M of two or more 
edges 

x(M; X, y) ^ cm{x + y) + higher terms in x, y , (1) 

the coefficients cm of x and y being the same positive integer. We first 
verified this on the computer for maps of up to ten edges. Since then 
Professor Tutte has pointed out to the author that (|l]) can be proved 
in general by recursively using Theorems 1 and 2 of Appendix A. 

The last case is equivalent (to within a re-definition of ^1,^2) to 
x{G) = 1 for all maps Q, i.e. to a Potts model with no interaction (the 
ideal lattice gas) . In this case the problem reduces to the enumeration 
of non-separable rooted planar maps. As the author discovered after 
solving this case, there is an extensive literature on this subject. |27, 
|29| , |33| , 36 , |2^, |8|, ^, ^ Explicit results have been obtained by Tutte 

12 1 They involve the 



(eq. 6.4 of Ref. [gg) and Brown and Tutte. [jTT^ 
intriguing hyper geometric function identities (|39 ), (^ 
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We also extend case 3 slightly, conjecturing the second-order cor- 
rections for X large when y = 1. 

In statistical mechanics we are particularly interested in large sys- 
tems, when we expect the bulk properties to be obtainable by taking 
an appropriate limit. Combinatorially this corresponds to the asymp- 
totic behaviour of the sum of dichromatic polynomials. In all our 
solved cases we find that the appropriate limit does indeed exist. 

For the regular planar lattices we know that the g-state Potts 
model has a first-order phase transition (i.e. a non-analyticity in the 
bulk properties) for q > qc, and a continuous transition for 1 < ^ < (?c; 
where qc = 4. |^, ^, 0, ^ It seems from numerical sudies that the 
random model of this paper behaves similarly, but with qc — 72. 

The outline of this paper is as follows. In the next section we 
define the dichromatic polynomial and the Potts model, and state their 
equivalence. Then we introduce the concept of non-separable rooted 
planar maps, and of the sum $ over their dichromatic polynomials. We 
give the functional relation satisfied by <I>. We then discuss the three 
explicitly solved cases, and go on to consider the asymptotic behaviour 
in the large- map limit. Finally we discuss the phase transition and 
numerical methods that can be used to study this. 

The dichromatic polynomial and the Potts 
model 

Let G be any connected graph (planar or not). Then the dichromatic 
polynomial is defined in as 

xiG; x,y)=Y.i^- 1)^(^)-^ {y - if(s)^E(s)-viS) (3) 
s 

where x, y are arbitrary variables, the sum is over all sub-graphs S of 
G, and C{S), E{S),V{S) are the numbers of connected components, 
edges and vertices, respectively, of 5. Since 5 and G share the same 
vertices, V{S) = V{G)- We have used Euler's relation: the number of 
circuits of 5 is 

P{S) = CiS) + E{S)-V{S) . (3) 

Thus for example the dichromatic polynomial of the graph of two 
vertices connected by a single edge is (x — 1) -|- (1) = x; for a triangle 
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of three vertices connected in pairs by three edges it is {x — l)'^ + 3{x — 
1) + 3(1) + (y- 1) = x2 + x + y. 

A g-state Potts model is defined on Q by associating with each 
vertex i a "spin" cTj that takes the values I,. . . ,q. Let K,v be two 
variables related by 

v = e^-l . (4) 
Then the partition function is 

ZiG;q,v) = ^exp{K^S{ai,aj)} , (5) 

where the outer sum is over all q^^^^ values of all the spins, the inner 
sum is over all edges of G, and 5{a, 6) = 1 if a = 6, 5{a, 6) = if 

The eqn. can be written as 

Z{g;q,v) = En{l + ^'^(^-'^j)} • (6) 

There are E{g) factors in the product, one for every edge of Q. Ex- 
panding the product, we obtain 2^*-^^ terms, each of which can be 
represented by a sub-graph 5 of ^ by including an edge in S if we 
take the v6{ai, aj)} term, leaving it out if we take the leading term 1. 

For each of the 2^^^^ terms thus obtained, we can readily perform 
the sum over the spins, giving 

=5:g^(^).^(^) , (7) 
s 

so we see that 

Z{g; q, v) = (x- 1) (y - l)^(^) x{0; x, y) (8) 

provided that 

q=ix-l)iy-l) , v = y-l . (9) 

Non-separable rooted planar maps 

A planar map (or simply a map) is a connected graph embedded in 
the surface of a sphere or closed plane, with no crossing edges. Not 
all graphs are planar, and those that are may correspond to more 
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than one planar map: for instance the maps in Figure ^ are different, 
but correspond to the same graph. If the vertices and faces are given 
distinctive labels, then a map can be specified uniquely by listing the 
anti-clockwise cyclic sequence of alternating vertices and faces (two of 
each) round each edge. Two maps are distinct if and only if their edge 
lists cannot be made the same by re-labelling the vertices and sites, 
and re-ordering the list. 




Figure 1: Two maps corresponding to the same graph. 



A map M is rooted if one edge A is chosen as the root-edge and 
given a direction. The vertex at the start of this directed edge is known 
as the root-vertex , the face on the left as the root-face. Two rooted 
maps are distinct if and only if their edge lists, root-edge, root-vertex 
and root-face cannot all be made the same by re-labelling the vertices 
and sites, and re-ordering the list. 

Finally, a graph is non-separable if there is no vertex such that the 
graph obtained by cutting it at this vertex (i.e. disconnecting from 
the vertex all its incident edges) is no longer connected. The non- 
separable rooted maps of up to three edges are shown in Figure ^. 
Their dichromatic polynomials, as defined by (^), are x,y,x + y, x^ ■\- 
x-\-y,x-\-y + y'^, respectively. The first two are the link-map and the 
loop-map: they are the smallest maps we shall consider. Note that we 
allow repeated edges, but only the loop-map contains a loop, and only 
the link-map has sites of valence one. 



Here we follow the original notation of Tutte[31| rather than Liu. 1 24] 
For a given map M, let m,n,i,j be the valency of the root-face, the 
valency of the root-vertex, the number of faces other than the root- 
face, and the number of vertices other than the root-vertex. Then 
the number of edges of M is i -|- j. We say that a map is "of type 
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Figure 2: Non-separable rooted maps with up to three edges. 
{m,n,i, j)" if it has these values of m,n,i,j. Define 

$ = $(ui,«2) = ^{ui,U2,Vi,V2,X,y) 

= J2<^2vi4x{M;x,y) , (10) 

M 

where the sum is over all distinct rooted non-separable maps M and 
m,n,i,j are the numbers defined above for each map M. Writing the 
contributions from the link and loop maps explicitly, this is 

^{ui,U2) = xulu2V2 + yuiulvi + ^'{ui,U2) , (11) 

where the sum $'(ni, U2) is further restricted to maps of two or more 
edges. Such maps also have i, j > 1 and m,n>2, so 

$'(txi,U2) = J2 N{m,n,i,j)uTu^v\vi , (12) 

m,n,i,j 

where 

l<m-l<i, l<n-l<i . (13) 

The coefHcient N{m, n, is the sum of the dichromatic polynomials 
over all distinct non-separable rooted maps of type {m, n, We also 
define coefficients N',P,R by 

$'(1,U2) = ^Ar'(n,z,j>^^;X2 , 
$'(1,1) = J2Pi''^)v\vi . (14) 

i,j 

00 

^'il,l,v,v,x,y) = Y.Rik)^" ■ 

k=2 
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Thus P{i,j) is the sum over all distinct non-separable rooted maps 
with i + 1 faces and j + 1 vertices; R{k) is the sum over all such maps 
with k edges. Clearly 

i+i 

N'{n,i,j) = ^ N{m,n,i,j) , 

m=2 

■i+l k-l 

P{i, j) = N'in, i, j) , R{k) = Y,P{i,k-i) (15) 

n=2 i=l 

From Figure § and the comments above, 
$'(tii,M2) = {x+y)u\u\viV2+{x^+x+y)u\u\viV2+{x+y+y'^)u\u\vlv2-\ 

all other terms being of order v\v2, with i + j > 3. 

Tutte shows that the concept of duality extends to rooted maps. 
If M* is the dual of M, and m*,n*,i* ,j* the corresponding values of 
m, n, i,j, then 

m=n,n='m,i=j,j=i , [lb) 

and 

x{M;x,y) = x{M*;y,x) . (17) 
The dual of a non-separable map is also non-separable, so 

^{ui,U2,vi,V2,x,y) = ^{u2,ui,V2,vi,y,x) . (18) 

We can use this duality to obtain the sum over n of N{m, n, from 
N'. 

Define 

^, X ni$(l,U2) - ^>(ui,U2) 

C{ui,U2) = , 

1 — Ul 

D{ui,U2) = ; . (19) 

1 - n2 



Then Liu[24| shows that $ satisfies the functional relation 

<^(U1,U2) = Xu\u2V2 + yUlU2Vl + 



.2, „, , „„ „2„, ^ UiU2VlC{ui,U2) 

1-C7(ni,l) 

ni^2^2-P(^l,^2) 

+ 1-I)(l,n2) • ^'"^ 
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For completeness we give the derivation of this result in Appendix A. 

We can solve this equation iteratively to obtain the contributions 
to $ from non-separable rooted maps of 1, 2, 3, ... edges. As we shall 
find in the next section, the numbers of such distinct maps are 2,1,2,6, 
22, 91, 408, 1938, 9614, 49335,.... (e.g. there are 91 non-separable 
rooted maps with six edges). 

As a check, we have enumerated all non-separable maps of up 
to 10 edges, and their various rootings and dichromatic polynomials, 
directly on the computer. We have verified that the results for ^ 
do in fact agree with those obtained from (pO[). We also find that 
the numbers of distinct unrooted maps are 2,1,2,3,6, 16, 42, 151, 596, 
2605,... The problem of counting unrooted maps has quite an extensive 
literature. Wormald [36|, and Liskovets and Walsh |22], established 
linear relations between the numbers of unrooted and rooted maps: 
both the above sequences of integers appear in Figure 1 of Ref. |22]. 
The relation between them is given in Appendix B. 



Case 1: x^y both large 

One case that can be easily handled is when x, y are both large and of 
the same order. In this case the coefficients in <^ of V1V2, for i,j>l 
and arbitrary ui,U2, are of order i + j + 5ii + 6ji — 2. To obtain them 
it is sufficient to set vi = wi/y,V2 = W2/X, take 11)1,102 to be fixed and 
of order one, and then to solve (|20| ) iteratively in inverse powers of x 
and y. To first order the the result is 

^/ N 2 2 ^ul:u%wiW2 _iU^ulwiW2 

<${ui,U2) = U1U2W2+U1U2W1 + y Vx 1 

1 — U1W2 1 — U2W1 

Continuing to second order, then setting ui = n2 = 1, we obtain 
<1>(1, 1) = wi + W2 + IJL1W1W2/X + wifj.2W2/y + fijwiwl/x'^ + 

ix\wip.2wl/{xy) + fiiwlfilw2/ {xy) + wlfii'W2/y'^ H (21) 

where for brevity we have written 1/(1 — Wi) as fii. From this and 
(0) it follows that if i = 1 or j = 1 

P{i,j) x^ii+y'5ji , 

ifi,j >2 

^'(i,j)~(i + J>W"'+rt^)(x^-5,2 + y^<J,2) , (22) 
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and 

i?(A;) ~ + . (23) 

Case 2: y small 

For any graph, the dichromatic polynomial vanishes when x = y = 
0. Consider the limit when x, y are small but non-zero, of the same 
order. Then any dichromatic polynomial is at most of this order. 
From and (ID, so are C and D. If we define Q = Q{ui,U2) = 



Qiui,U2,vi,V2) by 
^{ui,U2) = xulu2V2 + yuiulvi + {x + y)ululviV2Qiui,U2) , (24) 



then we find from ( |20D that x, y cancel out of the functional relation 
for Q, leaving 

X , , UiQ{ui,U2)-Q{l,U2) 

Q{Ui,U2) = 1 + U2V1 ; h 

Ui — 1 

U2Q{UI,U2)-Q{ui,l) 

U1V2 . (25) 

U2 - I 

This implies that, apart from the contributions of the link-map 
and the loop-map, $(tii, U2, wi, ^'2) is proportional to x + y when x 
and y are both small. As we state in the introduction, the stronger 
statement (|^) is true, for every non-separable map M of more than 
one edge. For all maps of two to five edges cm = 1, for six or seven 
edges it is 1 or 2, for eight to ten edges it is no bigger than the number 
of edges minus five (i.e. 3, 4, 5, respectively). The first map to have 
Cm > 1 is the tetrahedron (which is planar). 

The relation ([25| ) is a linear equation for Q{ui,U2)- If we regard 
Q{1,U2) and Q{ui,l) as known, we can solve it for Q{ui,U2)- The 
result is a rational expression with a denominator that is biquadratic 
in til, U2- If ^1, V2 are both small and U2 is of order one, the two zeros 
are when ui ~ 1 and when ui ~ l/{u2V2)- 

Suppose vi,V2 are both small, of order v. From (|l^ ), rn < i + j 
(for i + j > 2), so the coefficient of in the expansion of ^ is of 
order v"^ or smaller. It follows that the zero near ui = 1 must also be 
a zero of the numerator, since otherwise it would be a pole of <I> and 
the coefficient of would be of order unity. 

This gives a linear relation between Q{1,U2) and Q{ui,l) at this 
value of ui. We have not yet found a direct way to solve this relation. 
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but we have guessed explicit expressions for the coefficients in the 
expansion of Q{ui,U2)- In terms of the coefficients N,N',P of (|l2|), 
(0), they are: 

N{m, n,i,j) = (x + y){m — l)(n — l)A(m, n,i,j)x 
{i+j- 2)1 {i + j - m - + j - n - 1)! 



N'{n,i,j) 



{i - {i-n + 1)! {j - 1)! j! {j -m + 1)! ' 
{x + y){n - l)n(i + j - 2)! (i+j - 1)! {i + j - n - 1)! 



i\{i + l)\{i-n + l)\{j-2)\{j-l)\j 



p(- .^ ^ 2{x + y){i+j-2)\{i+j-l)\{i+3)\ 
^'^^ (.-l)!.!(. + l)!(j-l)!j!(i + l)! ' ^'^^ 

where 

A(m, n, i, j) = im + jn — mn + m + n — 2i — 2j . 



These formulae hold for all m,n,i,j satisfying (13) and i,j > 2. The 
last, for P{i,j), extends to i,j > 1, while the first two should be 
supplemented with 

N{m,n,i,l) = {x + y)5m2 5n,i+i , N{m,n,l,j) = {x + y)5n2Sm,j+i , 

N'{n,l,j) = {x + y)6n2 , N' {n,i,l) = {x + y)5n,i+i , 
for i,j > 1. 

It is straightforward to prove that these guesses are correct by di- 
rect substitution into (24) and (25). One needs the symmetry N{m, n, 
N{n,m, and the identity 

. .^ {i+j-m)B{m,n,i,j)N{m,n,i,j) 
}^N{k,n,i,j) = — — — , (27) 

where 

B(m,n,i,j) = (m — 2){m — l)i{i — n + l) + {j — 1)71(1771 — m — i+j + 1) . 



The identity (|27| ) can be verified by taking the difference between its 
m and m — 1 cases. 
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Case 3: enumeration of rooted non-separable 
maps 

If 

q = {x-l){y-l) = l . (28) 

Then the dependence on C{S) in (^) disappears and we immediately 
obtain (for all graphs) 

x(t?;x,y)=y^(^)/(2/-l)^(^)+i . (29) 

For our maps M, E{M) = i + j and V{M) = j + I, so 

x{M;x,y) =x^y' . (30) 

Hence x and y can be absorbed into V2 and vi and 

^'{ui,U2,vi,V2,x,y) = T{ui,U2,yvi,xv2) , (31) 

where 

T{UI,U2,VI,V2)=T{UI,U2)=Y.^TU244 ■ (32) 

M 

Thus the coefficient in T{ui,U2) of UiU2v\v2 simply the number of 
rooted non-separable maps of type {m,n,i,j). 

This case was solved originally by Brown and Tutte [|l^, and 
has been further investigated by a number of authors, including Ar- 
ques g, §, Bender |, |, Liskovets and Walsh |2|, Liu |3|, |2| 
and Wormald [^]. 

As a combinatoric problem this case is of interest in itself. From the 
point of view of statistical mechanics it is very simple, corresponding 
to the one-state Potts model or (as we show below) a q-state model 
with no interactions, i.e. the perfect gas. Even so, we hope to use it as 
the first step in examining the interacting Potts model of our random 
model, and we certainly do need to understand it if we are to use it 
as a starting point for general values of x and y. We therefore present 
the results here and outline the working in Appendix C. It involves 
some rather surprising identities. 

C and D simplify in the same way as <!>': C{ui,U2,vi,V2, x,y) = 
C{ui,U2,yvi,xv2) and D{ui,U2,vi,V2,x,y) = D{ui,U2,yvi,xv2). The 
definitions ([l9|) become 



A. s , UiT{l,U2,Vi,V2) -T{ui,U2,Vi,V2) 
C[Ui,U2,Vi,V2) = UiU2V2-\ 



1 — Ul 
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D{ui,U2,Vi,V2) = U1U2V1 + 



(33) 

U2T{ui, l,Vi,V2) - T{ui,U2,Vl,V2) 
I - U2 

Replacing 1^1,^2 in ( [20| ) by vi/y, V2/X and using these functions T, C, D, 
we find that x,y cancel out, except only that the last two terms (in- 
volving C and D, respectively) acquire factors l/y,l/x. This is the 
only place x and y now appear. Since the relation must hold for all 
x,y satisfying (p8|), and this can be written as 1/x + 1/y = 1, we 
obtain two relations from (pO|), namely 



UiU2ViC{ui,U2) UiU2V2D{ui,U2) 

^(^^1,^2) = — — = — — . (34) 

1-C;(ni,l) 1-D{1,U2) 

Other ways to reduce the dichromatic problem to simple enumera- 
tion are to allow x or y to tend to infinity, keeping the other fixed. For 
instance, when x is large the sum in (||) is dominated by sub-graphs 
S where C{S) attains its maximum value V{S). For all but the loop- 
map this means the sub-graph with no edges, so x(^i ^7 u) ~ £^nd 



we can absorb x into V2- We obtain (34), but without the last equality 
involving D. The resulting relation is still sufficient to iteratively de- 
termine T{ui,U2), and hence its desired power series coefficients, but 
it is not as full a description of the properties. 

These equations are fully solved in parametric form in Appendix 
C. In particular, we introduce two new variables p, r, defined in terms 
of f 1 , f 2 by 

vi=p{l-r)'^ , ^2 = ^(1-^)^ . (35) 



T{1,1) =pr{l-p-r) . (36) 



Then 



These equations are given in section 2 of ||T2[, our vi,V2, p,r being the 
X, y, u, V therein. As is noted in section 4 therein, products of powers 
of p and r (and hence of 1 — p and 1 — r) have remarkably simple 
expansions in powers of vi and V2- Let us introduce the standard 
notation 

{m)k = m(m + l)(m + 2) • • • (m + /c — 1) for k>0 , 

= 1 for A; = , (37) 
= l/[(m + A;)(m + A; + 1) • • • (m - 1)] for A; < , 
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where k is an integer. 

For positive integers j, define 

^, ,. 4(m + mj + mn)(2i + 2m + l),_i(2j + 2n + l)i_i 

iljl 

(38) 

This can be extended to i or j zero by taking an appropriate hmit, giv- 
ing P{m, n|0, j) = 2m(2m + (3{m, n\i, 0) = 2n(2n + 
and /3(m,n|0, 0) = 1. Thus P{m,n\i, j) is a polynomial in m and n, 
for all non-negative integers 

Then we find the following expansion in powers of vi , V2- 

CO CO 

pm^n ^ J2J2f3^„^^n\i,j)v\-^'^vi+^ . (39) 

i=0 j=0 

This elegant formula is true for all integers m, n, positive or negative, 
and can be extended to real and complex values of m, n. 
If vi = V2 = V, then r = p and 

oo 

p"" = Y^ aim, k)v^+'^ , (40) 

k=0 

where 

a{m,k) =2m{2k + 2m + l)''''^/k\ . (41) 

This is also an intriguing result. 

The coefficients P{i,j),R{k) in (|l§) are now easily determined: 

PH - ,.v 4(^ + 2j-2)!(2.+j-2)! 

(.-l)!(2z)!(,-l)!(2j)! ' ^"'^ 



and if X = y, 



2{3k - 3)! 



Thus P{i,j)/{y^'x^) is the number of non-separable rooted maps with 
i + 1 faces and j + 1 vertices, while R{k)/x'^ is the number with /c 
edges. These formulae were first given by Tutte and Brown p7| , |l2[. 

We have not found an explicit formula for the coefficients N' {n,i, j), 
N(m, n, i, j), but it does follow from Appendix C that they have a cer- 
tain structure. Let m' = m — l+6m2, n' = n — l+6n2, fn" = m— 5^2 — 2 
and n" = n — 6n2 — 2. Then 



N{m, n, i,j) 



Ay'x^ {2i - l)j_2m-n"+2 (2j - l)i-m"-2n+2 Gmn{hj) 

{i — n') \ [j — m')\ 
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^ ' (i-n + l)!(j-2)! ^ ^ 

Here Gmnihj) a-nd G'^{i,j) are polynomials in i,j of total degree 
2m" + 2n", 2n — 4, respectively. By "of total degree d", we mean 
that the sum of the powers of i and j in any term does not exceed d. 
In particular, G22{i,j) = 0, G23{i,j) = 1, G2,4(i, j) = 3(21 - 15i + 
2i2 - I7j + 5ij + 4j2), G^.^zihj) is of total degree 4, G2{i,j) = 1 and 
G'3(i, j) = 6z2 + 3i(5i - 11) + 12(i - - 2). 

The coefficients N(m, n, N'(n, are always finite. For small 
values of i,j the factors in ( ^4| ) other than Gmn{h j),G'ji{i, j) may 
become infinite, but in this case Gmnihj) or G'j^{i,j) vanishes and 
the coefficient can be evaluated by taking an appropriate limit. This 
is quite a strong constraint on the polynomials, but not apparently 
strong enough to determine them. 



The large-graph hmit 

In statistical mechanics we are usually interested in the "thermody- 
namic limit" , when the system becomes large in such a way that bulk 
properties, such as the total free energy or total entropy become pro- 
portional to the size of the system, and local or average properties, 
such as the density, tend to a limit. 

For any regular lattice the ratio of the number of faces to the 
number of vertices tends to a limit, determined by the valency. For 
the honeycomb, square and triangular lattices this ratio is ^, 1, 2, 
respectively. For any self-dual map it is one. This suggest that we 
should consider the limit when i, j become large, their ratio remaining 
fixed, non-zero and finite. This can be done in a symmetrical way by 
defining 

k = i + j , i = ka , j = k(3 , (45) 

and holding a, (3 fixed while A; — > oo. Note that a + (3 = 1. 

In all the cases discussed above we find that the coefficients 
N{m,n,i, j), N' {n,i, j), P{i,j), R{k) are asymptotically of the form 

constant x /C^/fc'^ , (46) 

i.e. their ratio to K,^ /k*^ tends to a non-zero limit as k ^ oo. This 
fC is the "partition function per site" and is in general a function of 
a, /?, X, y; c is an exponent, which for all our cases is a simple rational 
number. 
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Case 


Coeffs 


fC 


c 


1 


P 




-1 


1 


R 


max(x, y) 





2 


N, N', P 




9/2 


2 


R 


8 


4 


3 


N, N', P 




3 


3 


R 


27/4 


5/2 



Table 1: Values of /C and c for the coefficients N, N', P, R of the three solved cases. 
Define 

A = (l + a)2(l + /3)/(4a3) , 

^ = (i + «)(i + ;3)2/(4^3) _ (47) 

Then the results for IC and c are summarised in Table 1. 

For both cases 2 and 3, JC and c are the same for the coefficients 
N,N',P, provided we keep m,n fixed (or at least small compared 
with k) while we take the limit /c — > oo. Thus m, n affect only the pro- 
portionality factor in (^6[): from the viewpoint of statistical mechanics 
they are "boundary effects" corresponding to various weightings of the 
roofings. Indeed, it seems that the extra weights we have introduced 
by rooting the maps are also such a boundary effect: our result for R 
in case 3 can be compared with eqn. (3.21) of Zinn-Justin and Zuber, 
obtained by matrix model techniques. |3^] Allowing for an extra inte- 
gration, the two results agree to within a constant factor (independent 
of k). 

It is not actually necessary to calculate the coefficients to ob- 
tain /C. From ( [l^ and (^6|), the radius of convergence of the se- 
ries for <&'(1, 1, iJ, X, y) in powers of v is 1//C, so 1//C is the posi- 
tion of the closest singularity to the origin in the complex T;-plane of 
$'(1, 1, j;, y). For instance, for case 3 this is where the relation 
(|35|), with vi = V2, r = p, i.e. v = p{l — p)"^ ceases to be invertible. 
This is when dv/dp = 0, i.e. p = 1/3, v = 4/27, so /C = 27/4, as in 
Table 1. 

Similarly, if vi ^ V2 and p ^ r, then the singularity is when the 
Jacobian d{vi,V2)/d{p,r) vanishes, which is when 3pr+p + r = 1, i.e. 
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when r = (l — + 3p) . Thus 

' 2 i^(i+3p)2y \i+3p J ^ ^ 

For given a, /?, this is maximized when p is chosen so that 

p = a/{a + 2(3) , r = f3/{2a + (3) , (49) 

so then 

v\vi = \-'"^^-^f^ , (50) 

where X, fj, are defined as above. From Table 1, this is 1C~^, so the 
double series ( |l4|) for $'(1,1) then just fails to converge, as of course 
it should. 

Further discussion of the techniques available to handle the asymp- 
totic behaviour in case 3 is given by Bender and Richmond. P] 



A conjectured extension of case 3: y = 1 
and X large 

Ultimately it seems we should try to expand about case 3. We noted 
in that section that there are actually three sub-cases that reduce to 
the enumeration problem, namely x + y = xy, x large and y arbitrary, 
and y large and x arbitrary. We can consider an expansion about 
the second sub-case, expanding in inverse powers of x while keeping 
y fixed. If we choose y = 1, then we have observed a pattern in the 
coefficients P{i,j), namely that to second-next-to-leading order the 
RHS of (H) becomes 

,. 4(i + 2i-2)!(2z + j-2)! 
^ (i-l)!(2i)!(i-l)!(2j)! 

where 

^ ^ i{i - l)(2i - l)(i +i)(17 -9i + 3j - 3ij - 2f) 
{j + l)(2i + l)(2i + j - 2){2i+j - 3){2i+j - 4) 

We conjecture that this is correct for all i,j > l,2i + j > 4. It 
does have the expected behaviour ( |46| ) in the asymptotic limit. The 
partition function per site JC is modified, acquiring an extra factor 

a a + o? o?{oL + 2) 
^ X ^ ~ (3{a + 1)3 x2 + ■ ■ ■ 



l + ^ + ^^^^ + 0(.-3)) 

\ X 2x^ J 
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(remember that a + /3 = 1). To this order the exponent c remains 
at c = 3: there are no logarithmic terms in k that would imply a 
variation in this exponent. |18] 



Summary 



The functional relations developed by Tutte |31] and Liu |24] give 
equations that define the sum of dichromatic poynomials over rooted 
non-separable planar maps. This is the partition function of the Potts 
model, summed (i.e. randomized) over such maps. 

We are interested in the partial sum over maps with a given num- 
ber of sites and vertices (i + 1 and j + 1, respectively), or with a 
given number k = i + j of edges. We are particularly interested in 
the "thermodynamic" limit when i,j,k all become large, the ratios 
a = i/k,(3 = j/k remaining finite. We then expect the partial sums 



to behave asymptotically as ( 46 ) . The quantity /C therein is the "par- 



tition function per edge" . We expect it to depend on x, y, a, (3 (or just 
X, y if it is derived from R{k)), but not on the parameters m, n which 
relate to the way in which the maps are rooted. For this reason we 
would be happy to restrict attention to the case ui = ii2 = 1 (or any 
other fixed values of ui,U2), but we need to keep them arbitrary for 
the functional relations (^) as written to define the sum 

For the regular square lattice we know that the limiting function 
}C{x, y) has a singularity at x = y. For x > 3 (i.e. (7 > 4) this takes 
the form of a discontinuity in dfC/dx as one crosses the line x = y. 
For l<x<3(0<g<4) there is no discontinuity, but there is 
a non-analyticity. Physically this means that the Potts model has a 
phase transition at the self-dual point x = y, which is first order for 
(7 > 4, continuous for 1 < g < 4.[^, 37| 



The same is true for the triangular and honeycomb lattices, except 
that because these lattices are not self-dual, the transition is now at 
the point which maps to itself under a duality plus a star-triangle 
transformation. For the triangular (honeycomb) lattice this is when 
y'^ + y = x + 1 (x^ + X = y + 1).||] 

The obvious question is whether the /C of this paper has similar 
properties, or whether the transition is removed by the randomization 
introduced by summing over maps. For the regular lattices one can 
develop expansions for large x and y that make it quite clear that 
there is a first order transition. Unfortunately this technique fails for 
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our randomization. From Table 1 the exponent c is different for the 
three solved cases that we have presented. There are two explanations 
for this: either it varies continuously from case to case, or it varies 
discontinuously. We conjecture that the latter is true: all the values 
presented are integers or half-integers, and it seems likely that this is 
always true. Cases 1 and 2 are the large and small x, y limits, and it 
is likely that one is not allowed to interchange these limits with the 
large-/c asymptotic limit. This suggests that c has the values of case 
3 (namely 3 and 5/2) for all positive finite x and y, and this agrees 
with the second-order expansion of the previous section. 

So to investigate the phase transition problem we should focus on 
finite positive x, y and look to see if /C(x, y, a, /3) has non-analyticities 
for positive x, y, a, (3, in particular whether it is singular across the 
self-dual line x = y when a = p. 

In one sense the problem is solved: one can iteratively solve (20) in 
increasing powers of vi and V2, the computing needed growing poly- 
nomially with the number of terms. A.J. Guttmann and others have 
argued that this is effectively a solution: it is exponential growth 
(which is usual in statistical mechanical problems) that can make a 
problem intractable numerically. 

Of course the problem with this is that one needs to take the limit 
of high powers of vi and V2 in order to evaluate IC. Any finite powers 
will only give a numerical approximation to /C, and it is notoriously dif- 
ficult to determine non-analyticities from numerical approximations. 
One could focus on the case vi = V2 = v, thereby calculating R{k), 
and assign numerical values to x and y. This would still leave one 
with the problem of handling a series expansion in the three variables 
ui,U2,v, SO the calculation of R{k) would grow with k as k^. It may 
be possible to make progress in this way, but it would certainly be a 
major numerical exercise. 

One could reduce the problem to one of working with functions 
of a single variable ui or U2 (but not both) by generalizing the argu- 
ments given after (25), which would reduce the memory problems in 
a numerical calculation. It would of course be even better to remove 
both ui and U2, or even obtain explicit equations for the asymptotic 
behaviour in the limit k — >■ oo. This can be done in case 3: can it be 
done in general? 

Since writing these last paragraphs we have numerically investi- 
gated the sum of P{i,k — i) in ([l5|). If x = y, then P{i,j) = P{j,i), so 
the terms in the sum are symmetrical about i = k/2. For sufficiently 
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small X there is a single maximum, at i = A;/2, and this dominates 
the sum in the limit of k large. For sufficiently large x there are two 
maxima, symmetric about i = /c/2, which together dominate the sum. 
If X is slightly different from y, one will be larger than the other. As 
k increases, both will grow exponentially, and the larger will by itself 
dominate the sum. Hence one can define an order parameter: 



i=k~l 




This will be zero for x < x^, non-zero for x > Xc, where Xc is the value 
of X at which P{i, k — i) first becomes two-peaked. Numerical studies 
up to A; = 56 are reported in Appendix D. They indicate that Xc — 9.5, 
corresponding to Qc ~ 72. Above this value there is a discontinuity in 
the first derivative of the limiting value /C of Rik)^/^ across the line 
X = y in the (x, y) plane. This is a first-order transition. 

For the regular lattices, the singularity in the free energy of the self- 
dual Potts model at g = = 4 is essential (going like the exponential 
of —l/^/\q — qc\ ) '■ all derivatives exist and are continuous, but the 
function is non-analytic fs], . It would be interesting to determine if 
the same is true for this random model. 
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Appendix A 



Here we give the derivation of Liu's functional relation (2C). If yl is 
any edge of a graph Q that is not a loop or an isthmus, let Q'j^ be 
the graph obtained from Q by deleting the edge A, and t/^ the graph 
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obtained by identifying the two ends to form a single vertex. Then 
Tutte 1 31 1 gives the following theorem: 



THEOREM 1. 

x{G;x,y) = x(.GA,x,y) +x(.GA:X,y) . (Al) 

Also, if Q is the union of two subgraphs Tii and H2 having just 
one vertex in common (so G is separable), then Tutte also gives 

THEOREM 2. 

x{Q;x,y) = x(T-Ci;x,y)x(T-C2;x,y) . (A2) 

Tutte used these theorems to obtain a functional relation for the 
generating function of the sum of dichromatic polynomials over all 
rooted maps. Here we shall adapt his method to the sum over all 
non-separable rooted maps. 

The only non-separable maps with a loop or an isthmus are the 
first two graphs in Figure ^. Applying Theorem 1 to all other maps, 
we obtain 

$ = xu\u2V2 + yuiu^vi + Zi + Z2 , (A3) 

where 

Z, = J2x{M'A;x,y)uTu^v\vi , 

Z2 = Y.x{Ml;x,y)uTu^v{vi , (A4) 

the sums being over distinct non-separable rooted maps M with at 
least two edges. 

The map M'^ may be separable, but only in a very specific way: if 
s is the root vertex of M and t the other end of the edge A, there may 
be a sequence of nodes through which any route (after deletion of the 
edge A) from s to t must pass. These, and only these, are vertices at 
which Af^ can be separated. For instance, in Figure |^ all routes from 
s to t (not along A) must pass sequentially through the two nodes p 
and q. If / is the root face of M, and g the other face adjacent to A, 
then p and q are points where the faces / and g touch. 



20 




Figure 3: A map M that separates after removal of the root edge A. Each shaded 
region represents a non-separable map with an outer face of at least two edges. 



Consider the situation shown in Figure ^, where consist of 
three non-separable maps Mi, M2, M3 linked at p and q. We ask how 
many ways three such given rooted maps can be combined as in the 
Figure. For Mi we can take the root-vertex to be s, the root-face 
to be g, and the root-edge to be the first edge from s alongside g. 
Similarly, M2 has root-vertex p and root-face g; M3 has q and g. Let 
ma,na,ia, ja be the values of m,n,i,j for Mq,. Also, let the number 
of edges of Mq, adjacent to the face g be Tq,. Then nia > > 1, 
ria > 1. Also, m = mi + m2 + m-3 — ri — r2 — + 1, n = ni + 1, 
« = «i + ^2 + is + 1, j = ji + 32 + is- Using Theorem 2, the total 
contribution to Zi from the terms indicated by Figure |3| is therefore 

mi — 1 

UIU2VI J2 uT'-'^'u^^v\^vi'x{Mi-x,y) 

Ml ri=l 

r?i2 — 1 m-j — 1 

M2 r2=l M'i r3=l 

(A5) 

Note that U2 occurs only in the outer factor and the first sum. 

Strictly, the M^-sums should be restricted to non-separable maps 
Mq, with ttIq > 2. However, the only non-separable map with m < 2 
is the second map (the loop- map) in Figure with m = 1. If nia = 1, 
there are then no terms in the Tq summation, so the loop-map does 
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not contribute and we can extend the summations to all non-separable 
maps. 

Each summation over u^°'~^" gives a factor 

1 — Ui 

(Note that this vanishes for rria = 1-) Inserting these, using (|l^) and 
writing U2, fi, f2, y) simply as <I>(ni,-U2)! (|A5D reduces to 

UlU2VlC{ui,U2)C{ui,l)'^ , (A6) 



where C{ui,U2) is defined by (|19D. 

Figures like Figure ^, but with k + 1 separable pieces instead of 
three, give the same contribution ([A^), but with C(ni,l)^ replaced 
by C{ui, l)'^. Summing over A: = 0, 1, 2, . . ., we obtain 

Zi = uiU2ViC{ui,U2)/[l-C{ui,l)] . (A7) 

To calculate Z2 we note that Z2 is the dual of Zi and obtain 
immediately 

Z2 = UiU2V2D{ui,U2)/[l-D{l,U2)] , (A8) 



where D{ui,U2) is also defined by (|19|). 

From (lAST) , we therefore obtain the functional relation (l20|) . This is 



equation (4.17) of |24|, fi, u, x, y, z, t therein being our x, y, V2,vi,ui,U2- 



Appendix B 

Let R{k) be the number of rooted non-separable maps with k edges, 
and U{k) the number of unrooted ones. Then i?(l), -R(2), . . . = 2,1,2,6, 
22, 91, 408, 1938, 9614, 49335,... and U{1),U{2), . . . = 2,1,2,3,6, 16, 
42, 151, 596, 2605,... 

Also let <l){k) be the Euler totient function, i.e. the number of 
integers less than k that are relatively prime to k (including 1). Thus 
<^(l),(/.(2),... = 1,1,2,2,4,2,... 

Then Liskovets and Walsh ||2^ show that 

AkU{k) = 2R{k)+Y,Hk/j)(.9f -9j + 2)R{j) + kckR{k') , (Bl) 
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where 



Cfc = (A;+l), k' = {k + l)/2 if /c is odd , 
Cfc = (3A; - 4)/4, k' = k/2 if /c is even , (B2) 

the sum being over ah positive integers j that divide k and are less 
than k. 

Since this R{k) is known, being given by (^) with x = 1, the 
relation (|B1| ) defines U{k). The number of ways of rooting a given 
map is at most twice the number of edges (two directions for each 
edge), so 

R{k)<2kU{k) . 

Numerically we observe that 2kU{k)/R{k) decreases to unity expo- 
nentially fast, suggesting that this upper bound on R{k) is attained 
asymptotically. 



Appendix C 



To solve the equations (p4D, set M2 = 1 in the equation not involving 
D and write t{ui) for T{ui,l,vi,V2). We can eliminate the function 
C. Writing ui simply as n, we obtain the functional relation 



t{u) = [uvi + t{u)\ 



ut(l)-t(u) 

UV2 + 



(CI) 



1-u 

Writing t(l) as ti, this can be written as 

t{u)'^ + {1 — U + U^V2 — UV2 + UVl — Uti)t{u) + 

U'^ViV2 — 'U^V\V2 — U^Viti = . (C2) 

Given ti, this is a quadratic equation for t(u). Its discriminant is a 
quartic in u. We consider the situation when vi and V2 are both small 
and of order v. This quartic then has two zeros close to one, and two 
of order 1/v. 

The solution t{u) of the quadratic is certainly analytic for suffi- 
ciently small n, so has a convergent Taylor expansion. For a given map 
M (other than the link map) of type (m, n, i > 1 and j > m — 1, 
so i+j > m. The coefficient of u"* in the expansion of t{u) is therefore 
no bigger than order u™, so the radius of convergence is of order \/v. 



23 



It follows that the two zeros of the discriminant close to one must 
coincide, otherwise t{u) would have a square root singularity at them. 

This condition gives a fifth-degree equation for ti in terms of vi 
and V2- We do not write it down here as it is far better to proceed 
as follows. As a polynomial in u, the discriminant has the form 1 + 
■ • • + ^2^^. It must therefore be identical to the polynomial (1 — 
eti)^(l + 2fu + s^ti^), where e, /, s are parameters to be determined. 
Equating coefficients of u, . . . . ,u'^, we obtain four relations between 
vi,V2,ti, e, /, s, the last of which can be written as V2 = es. If we then 
define p, r so that e = 1 — p,s/e = r, we find / = r{2pr — p — 1) and 
that ^1,^2,^1 = T{1, 1) are given by (|35| ) and (^6|). 

We thus have a parametrisation in terms of the variables p, r. 
Given vi, f2, we can in principle solve for p, r and then obtain r(l, 1). 
We should choose the solution where p ~ vi,q ~ ?;2,T'(1,1) ~ viV2 
when vi,V2 are both small. Note that T(l, 1) is plainly a symmetric 
function of , f 2 , in agreement with duality. 

Writing u again as ui, we can go on to solve for t{ui) by replacing 
ui by wi, where 

.1 = . (C3) 

1 — r + r(l — p)wi 

As wi increases from to 1, ui also increases from to 1. 

The discriminant of the quadratic is now a perfect square and we 
find 

t{u,) = T{u, ^^^ = Prwl[l-p-r+pr'-pr{l- p)w,] ^^^^ 

1 — r + r(l — p)wi 

choosing the solution which vanishes when ui and wi are zero. 
If we also define W2 in terms of U2 by 

W2 {l-prw2) , . 

U2 = ■ — 7- ^ — , (C5 

1 — p + p[l — r)W2 

then from the duality symmetry 

prw'^[l-p-r+p'^r-pr{l-r)w2] 

T{l,U2) = ■ — 7- 7 . (C6) 

1 — p + p[l — r)W2 



It is now easy to solve ( |33D - ( p4| ) for the full function T{ui,U2) 
(the equations are linear in this function), giving 

T{ui,U2) = pr wlw2 {1 — prwi) {1 — prw2) x 
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(1 -p){l - r) - pr{l - p)wi -pr{l - r)w2 , . 

(1 — prwiW2) [1 — r + r(l — p)wi\ [1 — p + p{\ — r)w2] 

This completes the parametric solution of the functional relation 
for T{ui,U2), but it is still far from obvious that the coefficients in the 
expansion in powers of vi,V2 will be anything straightforward. In fact 
they are (at least for T{1, 1) ), as one can observe empirically. The 
key to the proof is the identities (|39|), (|40|). 

Let us look at the second (simpler) identity (^) first. The rhs is 
certainly some function of m: write it as P{m). If ( |35|) is true, then 
p is related to vi = V2 = v by v = p{l — p)'^ . Multiplying by p™, this 
implies the identity 

P{m + 3) - 2P{m + 2) + P{m + 1) - vP{m) = . (C8) 



This is easily verified from the series expansion ( ^OD of P{m). 

But ( p8| ) is a third-order difference equation. For integer m it 
follows that 

P{m) = Aip^ + A2P'^ + A^pf , (C9) 

where pi,P2,P3 are the three roots oi v = p{l —p)"^ and Ai,A2,A^ are 
some coefficients, independent of m. 

If we write (C9) down for three successive values of m, i.e. m,m + 



l,m + 2, we can solve the resulting three equations for 741,^42,^3 
(provided pi,P2,P3 are distinct, as they are for \v\ < 4/27). The 
result has the form 

Aj = [cijP{m) + C2jP{m + 1) + C3jP{m + 2)]/p™ , (CIO) 

where the coefficients cij, C2j, c^j are independent of m. 

For V sufficiently small, we can bound the coefficients in the ex- 
pansion (^) and show that P{rn) is of order f™. The three roots of 
the cubic are close to v, 1, 1; we choose pi = p v, p2 — Pz ~ I- 
Now let m ^ oo in ( plO ): the rhs vanishes exponentially for j = 2, 3, 



so A2 = A^ = 0. Hence P{m) is proportional to = p^. From 
(|40|), P(0) = 1, so P{m) = p"^ and we have proved (|40| ) for sufficientl 
small V. Both sides exist and are analytic for \v\ < 4/27, so it is true 
throughout this domain. 

When n = 1 the theorem implies that 

a result that can be obtained from equations (2.8.11) of Q. For all n, 
P(n) = t;"3F2(f ,^,^;n+i,n+l;^) , (C12) 
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3-^2 being the generalized hyper geometric function (section 9.14.1 of 
Ref. H). 

The first identity (|39|) can be proved in the same manner. One 
first ehminates r between the equation (^5|) to obtain ^1(1 — p)^ = 
p[{l — p)^ — ^2]^- This is fifth degree equation for p. When vi,V2 
are small the roots are approximately ui, 1,1, 1,1. Multiplying the 
equation by p"^ gives a fifth-order linear difference equation for the rhs 
P{m,n) of (pg|), as a function of m. We can verify that it is satisfied 
by the given series expansion. Proceeding as before, it follows that 
P{m,n)/p^ is independent of m. By symmetry, P{m,n)/{p'^r^) is 
independent of both m and n. We can verify that P(0, 0) = 1, so 
P{m,n) = p^r"" and we have proved ( |39| ) for sufficiently small vi,V2- 
Both sides are analytic within the radius of convergence of the double 
series, so the identity is true in that domain. 

Appendix D 

For X = y and k = 2h an even integer, let j{k) be the relative difference 
between the central value (with i = j = h) of P{i,j) and the adjacent 
value (with i = h - 1), so -f{k) = 1 - P{h - l,h + 1)/P{h,h). This 
is a measure of the curvature of the function P{i, k — i): positive for 
X < Xc, negative for x > Xc- For k finite we can define Xc to be the 
value of X for which j{k) vanishes. As k increases we expect it to tend 
to a limit, namely the bulk value of Xc referred to above. 

We have calculated Xc for = 6, 56 and give the results in Table 
2. They do indeed appear to be converging to a limit rather less than 
9.5, but it is hard to be more precise than this. One would expect such 
values to converge monotonically to a limit as an inverse non-integer 
power law in k, but these initially decrease, then increase to a maxi- 
mum at k = 48, and then start to decrease again. It looks as though 
even higher values of k are needed to examine the convergence. As a 
check of our numerical accuracy, we have performed the calculations 
to both 16-digit and 19-digit precision in Fortran (real*8 and real*12): 
the results agree to better than the accuracy of Table 2. 
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